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Abstract 

A systematic solution of a general model for Copper-Oxides reveals a line 
of transitions T = Tp(x) for x, the doping away from half- filling less than a 
critical value, to a phase with broken time reversal and rotational symmetry 
but with translational symmetry preserved. The single-particle spectrum in 
this phase is calculated to have a gap of dx2_y2 symmetry. The properties 
in this phase are compared to the properties of the so called "pseudo- 
gap phase" in the Copper-Oxides. The fluctuations towards this phase 
promote superconductivity which is either of the dx2_y2 or the generalised 
s-symmetry depending on whether the projected density of states at the 
chemical potential is larger in the dx2_y2 channel, as in the hole-doped 
copper-oxides, or in the s-wave channel. 
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INTRODUCTION 



A schematic phase-diagram of the Copper-Oxide (CuO) metals is shown in fig.(l). 
The superconducting region is surrounded by three regions with distinct properties: A 
region marked (III) with properties characterestic of a Fermi-hquid, a region marked (I) in 
which a Fermi-surface is discerned but in which the quasiparticle concept is inapphcable, 
and a region marked (II), the so-called pseudo-gap region, in which the concept of a 
Fermi-surface itself is lost. 

The topology of fig.(l) around the superconducting region is that expected around 
a quantum critical point in itinerant Fermions. The marginal Fermi-liquid (MFL) 
phenomenology with which many of the unusual properties in region (I) are understood 
assumes a scale-invariant low energy fluctuation spectrum characterestic of a quantum- 
critical point at X = Xc, the composition near the highest Tc. A region of Fermi-liquid 
is then expected for x > Xc at low temperatures, as in region (III), and a region with a 
broken symmetry for x < Xc at low temperatures. As in heavy-Fermion compounds, [Q 
a region of superconductivity is found at low temperatures peaked in the region around 
the quantum critical point. From this point of view, the crucial question in CuO metals 
is the symmetry of the phase in the region (II) below Tp(x). 

A systematic theory starting with a very general model for CuO provides an answer 
0] to this question. The region (II) in Fig. (1) is derived to be a phase in which a 
four-fold pattern of current flows in the ground state in each unit-cell as shown in Fig. 
(2). Four- fold symmetry around the copper-atoms as well as time- reversal symmetry are 
broken while their product as well as lattice translation symmetry are preserved. This 
phase will be referred to as the Circulating Current (CC) Phase. The properties of this 
phase were not studied in Ref. (4), which was mainly concerned with the fluctuations 
leading to the MFL properties in region (I) and the superconductivity induced by them. 

In the last few years, the properties in region (II) have become much clearer, thanks 
specially to Angle-Resolved Photoemission (ARPES) |^ and thermodynamic measure- 
ments. A pseudo-gap in the single particle spectra, consistent with dx2_y2 symmetry 
begins to develop at T < Tp(x) which is similar to the temperature at which the specific 
heat coefficient 7(T) |^ and the magnetic susceptibility x(T) |0] begin to decrease and 
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the resistivity and optical and Raman-spectra |T0| begin to drop below the MFL 
behavior of region (I). I calculate in this paper that the (CC) phase has a single-particle 
spectra with a gap of dx2_y2 symmetry, as observed, and of the right order of magnitude, 
and from which the other properties in region (II) follow. The principal features of the 
phase-diagram in fig.(l) then arise in systematic calculations from a single model. I 
suggest further experiments to confirm this identification. (Numerous ideas 0] already 
suggested for the pseudo-gap phase are discussed elsewhere fl^). 



THE MODEL 



The CC phase is a mean-field solution of a general Hamiltonian in the basis of three 



orbitals per unit cell, d, Px,Py: []T3 



H = K + h2 + Hg; 



K = ^ CdUdka + 2 tpdd]^^ (sx(k)pxk<7 + Sy(k)pyka) - 4tpp Sx(k)sy(k)p+]^^pyka + H.C. (2) 

k,(T 

Here a particular choice of the relative phases of the x and y orbitals in the unit cell has 
been made, Sx.y(k) = sin(kxa/2, kya/2) and for later, Cx,y(k) = cos((kxa/2, kya/2) and 
Sxy(k) = (sin^(kxa/2 + sin^(kya/2)). I consider the local interaction on the Cu and the O 
orbitals, 

Hint = X! ^dio- ndi-cr + Up (Upxicr T^pxi-a + Hpyio-npyi-o-) (3) 
i,(T 

and the nearest neighbor interaction between the Cu and the O orbitals, 

Hini = 2V ^ Cx(q) d^+^^ dk<x Pxk'_q^' Pxkv + X y (4) 

kk'q,cr<T' 

More general interactions do not change the essential results derived here. Throughout 
this paper (renormalised) energy difference ed between the Cu and the O orbitals is 
taken zero. It is important that in CuO, ed < O(tpd). Taking it as zero simplifies the 
calculations presented; the principal effect of a finite ed will be mentioned. 
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THE CIRCULATING CURRENT PHASE 



First, I derive some results of Ref. (4) in a simpler way by a simpler treatment of the 
large on-site repulsions in Eq. (3) and calculate the phase diagram of the CC phase. In 
the limit (Ud,Up) >> (tpd,tpp), a good mean-field approximation for low density of 
holes (or electrons) consists in replacing 

tpd ^ tpd = tp(j|x| 

tpp ^ tpp = tpp|x| (5) 

where |x| is the deviation from half-filling in the conduction band: x > for holes 
and X < for electrons. A more general treatment would consider separately the average 
occupation in the oxygen and copper orbitals and renormalize tpd, tpp accordingly. This is 
an unnecessary complication - no new physical principle is involved and small quantitative 
modifications are expected since the average occupations of Cu and O orbitals in CuO- 
metals are quite close (near optimum doping). 

A mean-field (non-superconducting) order parameter is sought which does not break 
translational symmetry (nor the rotational symmetry of spins). This means that the 
mean-field Hamiltonian is just a change of the coefficients in the kinetic energy opera- 
tor K. All the possible mean-field decompositions of the interactions except one only 
change the magnitude of the coefficients while preserving the symmetry. The only in- 
teresting mean-field decomposition comes from H-j^/ and yields the complex mean-field 
order parameter 

Rei^ ^ V/2 Y: Sx(k) (d+ p.k.) - s,(k) (d+ p,k.) (6) 

ko- 

The mean-field Hamiltonian itself is 



H^f = K - Re-'* J2 (sx(k) d+ p.k. - (k) d+ p,k.) + H.C. + ^ (7) 

k(T 

By a unitary transformation 

dko- — * dkcr, Pxko- Pxkcr e'^", Pykcr ^ Pykcr 6^^^ , (8) 

tan^x^y = ±Rsin /(2tpd ± Rcos 0) , 4tjd^ ^ = (2tpd ± Rcos0)^ + R^sin^c/), (9) 
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the phase difference is transferred completely to the px — Py bond, so that Hmf — > Hg, 
He = ^ 2 [dt„ (tpdxSx(k)pxka + tpdySy(k)pyk^)] - 4tpp e'^Sx(k)sy(k) p+ ^ pyk^ + H.C. + RV2V , 



(10) 



where 6 — 6y_ — 9y. It is apparent that only the phase difference of the two p — d bonds 
in the unit cell is gauge-invariant. With the most general interaction terms, the gauge- 
invariant phase is the change in phase going around the plaquette formed by a Cu-atom 
and any two of its oxygen neighbors, which differs from the above only in the value of 9. 

The energy of each of the three bands obtained by diagonalising Eq. (7) is changed 
by a finite (R, 0). But for any k the trace of the change in energy of the three bands is 
zero. The change in energy can therefore be expressed purely in terms of the change in 
energy 5 eck of the filled part of the (hole) conduction band. So the mean-field values Rq 
and 00 are determined by minimizing 

|^-2X:5eck(0,R)f(eck). (11) 

5eck has a complicated expression. To O(R^) and leading order in (tpp/tpd) it is 

R 2 ... 2.^ , 'l(tppS2(k)s2(k)) 1 



^eck^^l;^) I tpdS^y(l - l/4cosV) + ^ ' (1 " ^ ""^^ ^) I " (12) 

Prom Eq. (11) and (12) Rq 7^ at T = for 

2|x|tpd ^ ^ / 2 ,^ , 8tpps2(k)s2(k)\ 



< E + Y 

k<kF \ 



(13) 



^ k<kF V S2y(k) ) 

and 00 is vr/2 or — 7r/2. The symmetry of the transition is therefore of the Ising variety. 
For |x| << 1, Eq.(13) is satisfied only for x less than critical doping |xc|, 

|xc| ~ ^(V/(tpd)(0.25 + 0.5tpp/tpd). (14) 

|xc| defines the quantum critical points for bothe electron and hole dopings. An estimate 
of 9q = Ro/(2tpd) obtained by expanding Eq. (12) to 0(R**) is 

9l = O (tJd/2[V(tpp + tpd)])(xe - x)) . (15) 
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For ed 7^ 0, t^^ is replaced by tp^ + O(e^). With V, tpp,tpd, and ea of similar magnitude, 
as for CuO compounds, Xec,Xhc are therefore about 0.2. The line of transitions at finite 
temperatures varies with x as 

Tp/Ep- |x,-x|V2 ^^g^ 

This is to be identified with Tp(x) of fig.(l). The competition with antiferromagnetism 
near x = is not treated here. 

The energy of the (transverse) fluctuations in (p about 0o at long wavelengths is 
estimated from Eq. (12) to be 

^o = 0(eo'(tpp + tpd)/4). (17) 

The variation of fiq with q is estimated to be slow, of 0(qa)^. 

The eigenvectors of the states in the conduction band to leading order in 9 and tpp/tpd 

are 



|ckoa > +12^2^0^^ Sx(k)sy(k)|akoa > 
tpd 



where |ako > and |cko > are the nonbonding and conduction band states of k for = 
and Nk is the normalization factor. The term proportional to Oq leads to a difference in 
the phase between the Ipxko- > and |pyk(T> terms in Eq. (18) producing a time-reversal 
breaking state. Similarly expressions can be derived for the other two bands. The three 
together correspond to a current carrying state with pattern shown in Fig. (2). The 
wave function of states with 0o = — 7r/2 is (18) with — i^o in the second term. Note that 
states with 0o = 7r/2 and — 7r/2 are not orthogonal. But the two ground states formed 
from the Hartree-Fock product of such states are orthogonal. 

The CC phase breaks a rotational invariance of the Cu-0 lattice (besides the time- 
reversal invariance). So lattice defects, such as interstitials, dislocations or grain bound- 
aries couple as external fields do to a d = 2 Ising order parameter. [|15| No thermody- 
namically sharp transition is therefore possible at Tp(x). 
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ANISOTROPIC GAP IN THE NORMAL STATE 



In the mean-field approximation above, the energy of states in the conduction band 
shift in a k-dependent fashion, Eq. (12) with nothing special happening at the chemical 
potential. This is because the source of the instability considered is interband transitions. 
This solution will now be shown to be modified by the fluctuations which lead to scatter- 
ing among the conduction band states, (18), near the chemical potential. The fluctuations 
to be considered are transverse, i.e. of 0. whose energy Qq at long wave-lengths is given 
by Eq.(17). The Hamiltonian for such fluctuations, generated by operators (50q,50q, is 
quite generally obtained from the deviation of the starting Hamiltonian from Hmf . 



Hfluct = E + E g(k, k') c+,^ Ck,,, ((50k-k' + 50j _k) + H.C. (19) 

q kk'o- 

The coupling function to zeroeth order in tpp/tpd is calculated to be 

g(k,k') = 0otpd[Sx(k)Sx(k') - Sy(k)Sy(k')]/S.y(k). (20) 

For k k',this is ~ (cos(kxa) — cos(kya)). 

Consider now the renormalisation of the phase-fluctuation energy through Eq.(19). 
For small q 

= - 2nltl^9l{l - qVs) E ((cos(k.a) - cos(kya)/sxy)' /(ek+q - ek), (21) 

with the restriction on the sum that k < kf and k+q > kf. The self-energy of the fermions 
at k = kf and energy is ~ ln|ci; — figl suggesting an instability if the fluctuation energy 
at long wavelength approaches 0. The conclusion from Eq.(21) and Eq. (17) is that the 
phase fluctuations can be unstable over a substantial part of momentum-space. A way 
to restore stability is by a condensation of the phase fluctuations and a corresponding 
modification of the electronic energy. 

The condition for stability is derived by the ansatz that < S(j)q > and < C]^_^q g^Ckeo- > 
are finite for a range of q around zero [0. Eq. (19) then leads to new eigenstates which 
are annihilated/created by ■jka, ita- 

Ika = {Ckea + E V,k,k+qCk+qea)/Mk, (22) 



where is the normahsation. It will turn out that Uk,k+q 7^ only for states with 
|(e(k), e(k + q)) — /i| <^q ■ For such states momentum is not a good quantum number; 
I label the new states by k to indicate that the average momentum of such a state is k. 
In this sense, the problem retains translational invariance only on the average. 
From Eq. (19), 

1 ' 

< 50q >= -— ^ g(k, k + q) < C++q^e^Ck0<. > . (23) 
k 

I use the the Brillouin-Wigner (B-W) self- consistent approximation to get: 

< 4+ci,ea^k9a >= g(k, k + q) < 50q > /(Efc - Ek+q), (24) 

where E'j,s are the new one-particle eigenvalues to be determined. The B-W approxima- 
tion is exact in the limit that the number of states (k -|- q) coupled to a given state k is 
very large compared to 1, as here. 

Combining Eqs.(23) and (24) yields the self-consistency equations 

1 = ^ |g(k, k + q) I V(E.+, - E,), (25) 

k 

^k, k + q)^ g(k, k + q)/n; Yl g(k', k' - q)/(E,,_, - E',)^k', k' - q), (26) 

k'a' 

where 

$(fc, k+q) = {Ek+g - Ek) < c^+^0aCwa > ■ (27) 

In Eqs. (23)-(26), the sum is restricted to states k etc. such that e(k) is within about Qq 
of the chemical as required by the retarded nature of the Fermion-Boson interaction. 
The restriction on the sums is more clearly seen if Eq.(19) is used to generate an effective 
frequency dependent Fermion vertex 

g(k, k + q)g(k' - q, k')D^(q, uj)ci^^g^c^ea4'-q^sa'<^k'ea'- (28) 

k,k',q,(T,CT' 

where D<^(q, cu) is the propagator for the fluctuations which on the energy shell may be 
approximated by —l/fl'^ with a cut-off such that e(k'), e(k'-|-q) are both within about 
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of fi. A quadratic hamiltonian obtained by taking the expectation value of < c^^^^a'^kecr > 
etc. yields (26) using the B-W approximation. 

Eqs.(26),(27) have been obtained from the "direct" contraction of Eq (22). The 
exchange contraction yields a similar equation with a coefficient which, in the relevant 
channel, is (—1/4) of that on the right side of Eq (26). 

For an (approximate) solution of Eqs.(25),(26), note first that since q — > 0, g(k, k + 
q) ~ (cos(kxa) — cos(kya)), for $(k, k + q) ~ (cos(kxa) — cos(kya)) also. Second, 

E(k) ?a e(k) for |e(k) — /u| > Qq. A consistent solution for k near kf is then 

Ek = e], + D{k),k>kf 
Ek = ei,-D{k),k<kf 
D{k) = Do(cos(kxa) - cos(kya))^ (29) 

as may be verified by substitution. The value of Dq is evaluated to be of 0{9q tpjp(O)). 

The gap in the one particle spectrum at the chemical potential D(k) has a dx2„y2 
symmetry and a magnitude which can be estimated from Eqs. (15) to be x(xc— x)tp^(tpd+ 
tpp)/(Vtpp) which at say (xc — x)^ .05, and tpd ~ V ~ tpp ~ leV is ~ 20 meV. This has 
the right order of magnitude |^ . The differences of (29) from D-wave BCS single-particle 
spectrum are significant and have observable consequences discussed below. 

The leading decrease in energy due to the modification of the single-particle spectrum 
is ~ 6'q. Therefore, the nature of the transition remains unchanged,at least in mean-field 
theory; only quantitative changes are introduced in the conditions Eqs. (16) and (17) for 
the occurence of the CC phase. 



PROPERTIES IN THE CIRCULATING CURRENT PHASE 



The single particle density of states in the CC phase calculated from Eq. (29) is 



Pcc(c^) = p(0)-Arcsin 

71 





UJ 




UJ 


( 


D 




D 



< 1 



P(0), 



> 1 



(30) 



1 /2 

This increases as |ci;/D| ' for |u;/D| << 1 and unlike the d-wave superconductors, (which 
have a logarithmic singularity at = A, the superconducting gap), Pcc('^) is less than 



the "normal" density of states p(0) at all \u!\ < D. 

Eq. (30) should be compared with the single particle density of states measured 
in tunneling Jl^ and with the specific heat [H Cv and magnetic susceptibility x- 
The former shows, (just as Eq (30)), a diminuition in the single particle density of 
states for low energies at T < Tp(x) but show no rise above p(0) at finite energies until 
T < Tc(x), when the characteristic superconducting density of states appears. Cv in the 
CC phase is predicted ~ T^/^ and x ~ T^/^ for T << Tp(x). Due to the intervention of 
superconductivity, it is hard to test these power laws accurately. In the measured range 
I^J^ Tx/Cv is nearly independent of temperature as predicted here. x(T), measured 
more accurately than Cv(T), can be fit to T^/^. One can deduce the continuation to the 
T-dependence below Tc by invoking conservation of entropy on the Cv(T) measurements. 
Cv(T) ~ T^/^ for T << Tp(x) is then not inconsistent while Cv(T) ~ T clearly is. 

To obtain the spectral function A(k, u) measured by ARPES one needs besides 
Eq. (29), the self-energy of the single particle states (and the 'coherence factors'). The 
situation is in some ways similar to but in an important way different from the corre- 
sponding calculation for a d-wave superconductor (and completely different from 
s-wave superconductors). [|T3| In both cases, the bare polarizability xoi'i,^) (calculated 
from a single particle- hole bubble) is zero for uj < D(q), (A(q) for the superconductors). 
The lowest energy single particle scattering for momentum q occurs by an intermediate 
one-particle state near the zero of the gap. Therefore the threshold for single-particle 
scattering of a state at q is also D(q). For u > D(q), Xoi'i) with E^'s given by Eq. (29) is 
proportional to j^^. So the renormalized xi'i,^) for ^ > D(q) is expected to be similar 
to the normal state above Tp(x), i.e. of the marginal Fermi-liquid form. The single 
particle self-energy needs to be calculated in detail but its general form is evident: For 
(u;, T) << D(q) it is exponentially small, but for {uj,T) > D(q) it returns to the value 
Sn(co',q, T) without the pseudo-gap. So 



using Eq. (29) and (31) and the marginal Fermi-liquid form for Sn is plotted in fig. (3) 




(31) 



The spectral function at the Fermi- wave- vectors kf defined by E(kf) = ±D(kf) calculated 



for a few temperatures [|20|. A pseudogap in the direction kf appears below T ^ D(kf) 
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producing the illusion of 'Fermi-arcs' shrinking as temperature decreases. The lineshape 
in fig. (3) at low energies, in the pseudo-gap region, fits the experimental curves |Q 
well within the experimental resolution. A prediction following from the results of the 
previous section is that D ~ (xc — x) and also ~ (Tp(x) — T). 

The single particle spectrum in d = 2 d-wave superconductor |]rB[ only differs from 
the above because of a singularity in Xo(Q) ^) uj 2 A and q 2kf in the directions 
connecting the maximums of the gap. This leads to a pole in ReS(ti;, q) and therefore a 
sharp peak in A(c(j, q) at q on the Fermi-surface close to the maximum of the gap followed 
by a continuum starting at A. No such sharp feature is to be expected in the pseudo-gap 
region because no singularity in xoi.'^i^) exists. 

SUPERCONDUCTIVITY 

I discuss here the superconducting instability in the quantum fluctuation regime (Re- 
gion I in Fig. (1)) in which 9q = {] and effective electron-electron interaction is through 
fluctuation in 9. (The calculation in the CC phase is similar with the added complication 
that one must do pairing calculations in a gapped electronic structure.) The coupling of 
the 9 fluctuations to the Fermions has the same functional form as g(k, k') and leads to 
an effective pairing interaction in the spin-singlet channel: 

Hpair = E I V(k, kf D,(k - k', ^)c+ ,clk',-<.c_k,-aCk,. (32) 
k,k' 

Here De is the propagator of the 9 fluctuations, which is essentially momentum indepen- 
dent and has a frequency dependence of a scale-invariant form |^ . We may approximate it 
below its cut-off frequency by — (cUc)"^ The coupling function in Eq (32) is then expanded 
in a separable form in lattice-harmonics as follows: 

Vpair(k,k') = -VV(4cUe) Egi^i(k)^i(k') (33) 

i 

where i are the generalised s, dx2_y2, and dxy channels: 

r;,(k) = s2^(k),gs = l/2 
?7x2-y2(k) = l/2[cos(kxa) - cos(kya)], gx2_y2 = 1/2 

r7xy(k) = Sx(k)sy(k),gxy = -2. (34) 

11 



The methods of Ref.(3), may be used now to deduce that pairing is possible only 
generalized s-wave and the dx2_y2 channels. The choice between the channels is deter- 
mined by the band-structure near the Fermi- surface. In the hole-doped cuprates dx2_y2 
is favored because in this symmetry the gap occurs over regions of Fermi-surface where 
the bands-disperse the least (because of the proximity to the Van Hove critical point) 
thereby reducing the energy optimally. If fi is far away from the van Hove singularity, 
as in electron-doped cuprates the projected generalised s-wave density of states is 
expected to be larger than the projected dx2_y2 density of states and s-wave pairing is 
favored. 



CONCLUDING REMARKS 



Although, I have presented a systematic theory in agreement with the principal ex- 
perimental results, one can be confident of the applicability of the theory only if the 
structure in Fig. (2) is observed. I have suggested |^ elastic polarized neutron scattering 
to observe the orbitals moments in Fig. (2). Another way to test Fig. (2) is through the 
polarization dependence of ARPES in single-domain samples. Some other predictions 
are mentioned in the text. 

Improvements could be made in the theory, especially the derivation of the gap in the 
Circulating-current phase. Note that the conditions (26,27) are equivalent to requiring 
that the phase-fluctuation modes have zero-energy as q — > 0, which appears to give an 
xy-symmetry character to the circulating current state. If true, no specific heat anomaly 
is found even in the pure limit (except for three-dimensional couplings). The ansatz, 
Eq.(22), by which the problem has been diagonalised may be looked upon as admixing 
states with 0o = — vr/2 in to states with (po = tt/2. It is possible that a better treatment 
would express the one-paticle states as topological excitation of fermions bound to co- 
moving current fluctuation. These matters require further investiagation. 

I wish to thank E. Abrahams, B. Batlogg, P. Majumdar, and A. Sengupta for useful 
discussions. 
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FIGURE CAPTIONS 



Fig. (1): Generic phase-diagram of the cuprates for hole-doping. Not shown is a low tem- 
perature "insulating phase" in region II due to disorder. 

Fig. (2): Current pattern predicted in phase II of Fig. (1). 

Fig. (3): The spectral function at the Fermi- wave- vectors as a function of energy below the 
chemical potential normalized to the gap in the direction of that wave- vector D(kf) 
for temperatures T = nD(kf) / (An) for n = 1, 2, 3, 4, 5. 
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